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Abstract 

The main physical result of this paper are exact analytical solu- 
tions of the heavenly equation, of importance in the general theory of 
relativity. These solutions are not invariant under any subgroup of the 
symmetry group of the equation. The main mathematical result is a 
new method of obtaining noninvariant solutions of partial differential 
equations with infinite dimensional symmetry groups. The method 
involves the compatibility of the given equations with a differential 
constraint, which is automorphic under a specific symmetry subgroup, 
the latter acting transitively on the submanifold of the common so- 
lutions. By studying the integrability of the resulting conditions, one 
can provide an explicit foliation of the entire solution manifold of the 
considered equations. 
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1 Introduction 



An important problem for partial differential equations invariant with re- 
spect to an infinite Lie group is to obtain non-invariant solutions that admit 
no continuous symmetries of the equations. In our opinion, the old approach 
of S. Lie 0] developed by Vessiot and in modern form by Ovsiannikov p|, 
which we call group foliation, is an adequate tool for treating this problem 
in the framework of Lie theory. According to this method we foliate the 
solution space of the equations in question into orbits, choosing for the foli- 
ation an infinite-dimensional symmetry group. Each orbit is determined by 
the automorphic system joined to the original equations and considered as 
invariant differential constraints. Due to the automorphic property of this 
system, any of its solutions can be obtained from any other solution by a 
transformation of the chosen symmetry subgroup. This symmetry property 
makes the automorphic system completely integrable if only one of its solu- 
tions can be obtained. The collection of orbits of all solutions of the original 
equations is determined by the resolving system. Thus the problem reduces 
to obtaining as many particular solutions of the resolving system as possible. 
Each of them will fix a particular automorphic system and the corresponding 
orbit in the solution space of original equations. 

Group theory is usually used to obtain invariant solutions. Here we show 
that it also provides a mechanism for obtaining non-invariant solutions. We 
give examples of such solutions as an application of the method. 

In this paper we further develop the method of group foliation by in- 
troducing a procedure of invariant integration. It is used for reconstructing 
the solution of the original equation corresponding to the known particu- 
lar solution of the resolving system. We apply the method for obtaining 
non-invariant solutions of the 'heavenly" equation 

UXX + Uyy = K{e'')u (1.1) 

where k = ±1 and the unknown u depends on the time t and two space 
variables x and y. Here and further subscripts of u denote partial differen- 
tiation with respect to corresponding variables. This equation formally is a 
continuous version of the Toda lattice . It appears in various physical the- 
ories, like the theory of area preserving diffeomorphisms , in the theory of 
the so-called gravitational instantons and in the general theory of relativ- 
ity 1]^. In this context it describes self-dual Einstein spaces with Euclidean 
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signature with one rotational Killing vector. Moreover it is a completely 
integrable system in the sense of the existence of a Lax pair 

The outline of the method is the following. We determine the total group 
of point symmetries of the heavenly equation. For the group foliation we 
choose its infinite subgroup of conformal transformations. We compute dif- 
ferential invariants of this subgroup up to the second order inclusively and 
obtain 5 functionally independent differential invariants. On account of the 
heavenly equation we are left with 4 invariants. We choose three of them 
as new independent variables, the same number as in the heavenly equation, 
and one is left for the new unknown. 

We obtain three first order operators of invariant differentiation defined 
by the property that acting on a differential invariant they produce again a 
differential invariant. These operators are determined by the condition that 
they should commute with an arbitrary prolongation of any element of the 
infinite symmetry Lie algebra chosen for the foliation. 

Extensive use of operators of invariant differentiation and their commu- 
tator algebra for formulating the resolving system is a new feature of the 
method suggested by one of the authors (M.B.S.) in a recent article on the 
complex Monge- Ampere equation We derive the resolving system as a 
set of compatibility conditions for the heavenly equation and its automorphic 
system, using invariant cross-differentiation. Then we formulate the resolving 
system in terms of the commutator algebra of operators of invariant differ- 
entiation by discovering the fact that this algebra together with its Jacobi 
identities, projected on the solution manifold of the considered equation in 
the space of differential invariants, is equivalent to the resolving system. 

We show how an Ansatz simplifying the commutator algebra of opera- 
tors of invariant differentiation leads to a particular class of solutions of the 
resolving system. Then we use invariant integration to obtain the corre- 
sponding solution of the heavenly equation and prove that this solution is 
non-invariant. 
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2 Lie group of point symmetries and 
differential invariants 



It is convenient to work with the heavenly equation using the complex coor- 
dinates z = {x + iy)/2, z = [x — iy)/2 

u.-z = f^ientf (2.1) 

A standard calculation of the total symmetry group of the heavenly equation 
gives the following result for the symmetry generators of all one-parameter 



subgroups Wl 



T = dt, G = tdt + 29„, 
X, = a{z)d, + a{z)d, - {a\z) + a\z))d^, (2.2) 

where T is the generator of translations in t, G is the generator of a dilation 
of time accompanied by a shift of u: t = ie'^ , u = u+2t and Xa is a generator 
of the conformal transformations 

z = (l){z), z = (t){z), u{z,z,t) = u{z,z,t) — \n{(f)\z)4>'{z)), (2.3) 

where a{z) and (j){z) are arbitrary holomorphic functions of z (see also 0]). 

The Lie algebra of the symmetry generators is determined by the com- 
mutation relations 

[T, G] = T, [T, Xa] = 0, [G, Xa] = 0, [Xa, Xt] = Xat'-ba', (2.4) 

which show that the generators Xa of conformal transformations form an 
infinite-dimensional subalgebra. 

We choose for the group foliation the corresponding infinite symmetry 
subgroup of all holomorphic transformations in z, i.e. the conformal group. 
Differential invariants of this group are the invariants of all its generators Xa 
of the form ( p.2|) in the prolongation spaces. This means that they can depend 
on independent variables, the unknowns and also on the partial derivatives 
of the unknowns allowed by the order of the prolongation. The order N 
of the differential invariant is defined as the order of the highest derivative 
which this invariant depends on. The determining equation for differential 
invariants $ of the order N < 2 has the form 

1($) = 0, (2.5) 
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where X is the second prolongation of the generator Xa (|2.2| ) of the conformal 
group defined by the standard prolongation formulae 

2 

X = adz + o.dz — [a' + a!) (9„ — [a" + a'u^) du^ — {a" + a'uz) 

- {a" + a"uz + 2a!uzz) du,, - {a" + al'ui + ^a'u^z) du,, 

- a'uztduzt - a'uztdu-zt - {a' + a') Uzzdu,,, (2.6) 

where a = a{z) and a = a{z). 

The integration of eg . (12 .51 ) gives 5 functionally independent differential 
invariants up to the second order inclusively 

t, ut, utt, p = e'^'Uzz, V = e'^'uztu-zt (2.7) 

and all of them turn out to be real. This allows us to express the heavenly 
equation ( |2.1|) solely in terms of the differential invariants 

Utt = up- uf. (2.8) 



3 Operators of invariant differentiation and a 
basis of differential invariants 

Operators of invariant differentiation are linear combinations of total deriva- 
tive operators with respect to independent variables. Their coefficients de- 
pend on local coordinates of the prolongation space. They are defined by 
the special property that, acting on any (differential) invariant, they map 
it again into a differential invariant. Being first-order differential operators, 
they raise the order of a differential invariant by one. Invariance requires 
that these differential operators commute with any infinitely prolonged gen- 
erator Xa ( p.2|) of the conformal symmetry group. It is obvious (see par. 
24.2 for a complete proof) that the total number of independent operators of 
invariant differentiation is equal to the number of total derivative operators, 
that is to the number of independent variables (which is three in the present 
case). 

We look for operators of invariant differentiation in the form 

3 

6 = \,Dt + \2Dz + X^D-z = J2 (3.1) 

1=1 
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where Di = Df, D2 = Dz, D3 = Dg are operators of total derivatives with 

respect to the subscripts. We look for the coefficients Aj satisfying the condi- 

00 

tion of commutativity of 6 with the infinite prolongation X of the generator 
Xa ( |2.2| ). It can be decomposed as the sum of the infinite prolongation of the 
symmetry generator in the evolution form X and the linear combination 
of the total derivative operators 



X = X + Y, ^^Dj = X + a{z)D, + a{z)D-,, 
where from the form of Xa we take 



(3.2) 



(3.3) 



The generator X in the evolution form commutes with all the total derivatives 



Di,X =0 and hence we have the standard commutation relation 



(3.4) 



We use it in the determining relation for operators of invariant differentiation 
6,X] =E[a.A,x] =E(E^^A[e]/^,-X(A,)A,) 

i i j 

= E(E^.^.r]-^(A.) A = 0. (3.5) 



The final equation for the coefficients Aj of the operators of invariant differ- 
entiation (see eq. (24.2.3) of P) is 



x{\) = Y.\^Ai' 



(3.6) 



Using ( p.3| ) and restricting ourselves to the second prolongation X of the 
symmetry generator, the equation ( |3.6| ) leads to 



l(Ai) = 0, 1(A2) = \2a{z), 1(A3) = A3a'(^) 



(3.7) 
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where primes denote derivatives and 

Aj U, tij, U^, U^i U^ti U^ti U^zi '^^zzy '^zz)- 

2 

Here X is the second prolongation of the generator Xa of the conformal group 
defined by eq. (|2.6| ). 

Equations ( |3.7D are easily solved by the method of characteristics and we 
choose 3 simplest linearly independent solutions for the coefficients Aj of the 
three operators of invariant differentiation 

A} = 1, A^ = 0, A^ = 0, \l = 0, \l = e-'^Uzt, Xj = 0, 

A? = 0, Ai = 0, Ai = e-^uzf (3.8) 

From here we obtain a basis for the operators of invariant differentiation 

S = Dt, A = e-''u-ztDz, A = e-''UztD-z. (3.9) 

The basis of differential invariants is defined as a minimal finite set of 
invariants of a symmetry group from which any other differential invariant 
of this group can be obtained by a finite number of invariant differentiations 
and operations of taking composite functions. The proof of the existence and 
finiteness of the basis was given by Tresse and in a more modern form 
by Ovsiannikov 0. 

In our example the basis of differential invariants is formed by the set of 
three invariants t,Ut,p, while two other invariants Uu and rj of eq.( p.7|) are 
given by the relations 

utt = S{ut), f] = e~'^UztUzt = A{ut) = A{ut). (3.10) 

All other functionally independent higher-order invariants can be obtained 
by acting with operators of invariant differentiation on the basis {t, Ut, p}. In 
particular, the following third-order invariants generated from the 2nd-order 
invariant p by invariant differentiations will be involved in our construction 

a = A(p), a = A(p), r = 6{p) = pt- (3.11) 

The operators of invariant differentiation form the commutator algebra 

[5, A]=(^K^- 3u}j A, [6, A]=(^K^- 3ut^ A 
[A, A] = - {u,p + r)\A- - {utp + t)\ a (3.12) 
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which form a Lie algebra over the field of invariants of the conformal group, 
in agreement with Ovsiannikov's lemma 24.2 [Q. 

The commutator algebra is simplified by introducing two new operators of 
invariant differentiation Y and Y instead of A and A and two new variables 
A and A instead of a and a, defined by 

A = r/F, A = r/F, a = rjX, a = r]X, (3.13) 

and becomes 

[5, Y] = I^kX - 3u, - ^ j Y, [5, Y] = (^kX - 3n, - ^ j F, 

[Y,Y]=^-^^^P±^(Y-Y). (3.14) 

7] ^ / 



Equations (|3.10|) and (|3.11| ) imply the following properties of the operators 
Y and Y 

Y{ut) = Y{ut) = l, Y{p)=X, Y{p) = X. (3.15) 



4 Automorphic and resolving equations 

We have four independent differential invariants t,Ut,p,T] on the solution 
manifold of the heavenly equation ( p.8|) . We choose three of them t,Ut,p 
as new invariant independent variables, the same number as in the original 



equation ( p.l|) , and consider the fourth one ?7 as a real function F of these 
three 

7] = F{t,ut,p) UztUzte''' = F (t,ut,Uzze~'') , (4.1) 

which gives us the general form of the automorphic equation, i.e. invariant 
differential constraint. 

Our next task is to derive the resolving equations for the heavenly equa- 
tion. This will account for all integrability conditions of the system ( p.8|) , 
( |4.1| ) in an explicitly invariant form. If we pick a particular solution of this 
resolving system for F and use it in the right-hand side of ( [4.1|) , then the lat- 
ter equation will possess the automorphic property: each solution of it can 
be obtained from any other solution by an appropriate conformal symmetry 
transformation. 



We consider the automorphic equation ( [4.1| ) divided by F in the form 

Y{ut) = 1, (4.2) 
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and the heavenly equation ( |2.8| ) in the form 



5{ut) = Kp- uf. (4.3) 



We put ?7 = F in the definitions ( p.l3| ) of Y and Y and in their commutation 



relations ( |3.14 ). The integrability condition for the system ( [4 .21) and (|4.3| ) is 
obtained by the invariant cross- differentiation with 6 and Y with the use of 
their commutation relation ( p.l4| ) 

S{F) = [k{X + A) - 5ut]F. (4.4) 

Since this equation involves A and A we use their definitions in eq.( |3.15| ) 

Yip) = A, Yip) = A (4.5) 

and obtain the integrability condition for these two equations by the invari- 
ant cross-differentiation by Y and Y using their commutation relation from 
eq-(ra 

F(F(A) - F(A)) = iutp + r)(A - A). (4.6) 
This equation contains r, so we use its definition ( p.ll| ) 

5ip) = T. (4.7) 

Using the invariant cross-differentiation with Y or Y and 6, we obtain the 
compatibility conditions of eq.( [4.7| ) with each of equations ( [4.5| ) 



5(A) = r(r) + 2utX - kX^ (4.8) 

and 

5(A) = r(r) + 2utX - kX^. (4.9) 

These are complex conjugate to each other. There is one more differential 
consequence of the obtained resolving equations. This is the integrability 
condition of the equation ( [4.6|) solved with respect to V(A) together with the 
equation ( [4.9|) . It is obtained by the invariant cross-differentiation of these 
equations by 5 and Y . Using the other resolving equations it can be brought 
to the form 

F(r(A)+F(A)) = -(^itp + r)(A + A) 
+ 2/t[(5(r) +2F + AutT + Kp^ + 2u^^p] . (4. 10) 
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The resolving equations ( |4.4| ), ([4.6|), (|4.8|), ( [4.9| ) and ( |4.1(J| ) form a closed 
resolving system if we assume that not only the 2nd-order differential in- 
variant 7] = F, but also the 3rd-order differential invariants A, A and r are 
functions of t, Ut, p. They should be regarded as additional unknowns in these 
equations, so the resolving system consists of 5 partial differential equations 
with 4 unknowns F, A, A and r. The operators of invariant differentiation 
are projected on the solution manifold of the heavenly equation and on the 
space of differential invariants treated as new independent variables. We keep 
the same notation for the projected operators of invariant differentiation and 
write them in the form 
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dt + {np - ul)dut + Tdp, 



Y = dut + Xdp, 



(4.11) 



Here we have used the following properties of these operators 



6{t) = l, 6{ut) = Kp-ul 6{p) = T (4.12) 
Y{t)=Y{t) = 0, Y{ut) = Y{ut) = l, Y{p) = X, Y{p) = X, 



which follow from their definitions, equations ( |3.10| ), ( |3.11| ), ( |3.15| ) and the 
heavenly equation in the form ([4.3|). If we used for the operators of invariant 



differentiation 5,Y,Y the formulae ( [4.11 ) in the resolving equations 
( |4.6| ), ( |4.8| ), (|49|) and ( [4.10|) , then we would obtain the resolving system in 
an explicit form as a system of 5 first-order PDEs with 4 unknowns F, A, A, r 
and 3 independent variables t,Ut,p. This system is passive, i.e. it has no 
further algebraically independent first-order integrability conditions. 



The commutator relations ( 3.14 ) were satisfied identically by the opera- 
tors of invariant differentiation. On the contrary, for the projected operators 
( |4.11| ) these commutation relations and even the Jacobi identity 



[S, [Y,Y]] + [Y,[Y,5]] + [Y, [S,Y]]=0 



(4.13) 



are not identically satisfied, but only on account of the resolving equations. 
It is easy to check that even a stronger statement is valid. 

Theorem 1 The commutator algebra ( p.l4| ) of the operators of invariant 
differentiation S,Y,Y, together with the Jacobi identity ( [4.13| ), is equivalent 
to the resolving system for the heavenly equation and hence provides a com- 
mutator representation for this system. 
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This theorem means that the complete set of the resolving equations is en- 
coded in the commutator algebra of the operators of invariant differentiation 
and provides the easiest way to derive the resolving system. In Section 6 we 
shall show how the commutator representation of the resolving system can 
lead to a useful Ansatz for solving this system. 



5 Invariant and non-invariant solutions 

Invariant solutions are defined as solutions that are invariant with respect 
to a symmetry subgroup of the equation. Non-invariant solutions are those 
solutions which are not invariant with respect to any one-parameter sym- 
metry group of the equation. We present here a simple derivation of the 
infinitesimal criterion of invariance of solutions. 

Consider a general form of the generator of a one-parameter symmetry of 
the heavenly equation as a linear combination of symmetry generators ( |2.2|) 
with arbitrary real constant coefficients a and (3 

X = adt + P {tdt + 2du) + a{z)d, + a{z)d, - {a\z) + a\z)) 9„ (5.1) 

where a{z) is an arbitrary holomorphic function. The infinitesimal criterion 
for the invariance of the solution u = f{z, z,t) with respect to the generator 
X has the general form (see par. 19.2.1 of [^) 

X(/-n)|„=; = 0, (5.2) 

which for X defined by the formula ( |5.lD becomes 

(a + (3t)ft + a{z)f, + a{z)h = 2/3- a'{z) - a'iz). (5.3) 

The invariance criterion can be summed up as follows. 

Proposition 1 // there exists a holomorphic function a{z) and constants a 
and /3, not all equal to zero, such that the equation ( \5.3i ) is satisfied, then the 
solution u = f{z,z,t) is invariant. Otherwise this solution is non-invariant. 

From this proposition one can derive some criteria for the non-invariance 
of solutions. For example, we consider the case when a = and (3 = so that 
equation (|5.3|) is a criterion of conformal invariance. The general solution of 
eq.([5.3|) in this case has the form 



u = In f{^,t) -\na{z) -\na{z) (5.4) 
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where 

The invariant p defined by eq.( p.7|) becomes 

P 



P = ^^^^ (5.6) 



and the invariants a and a, defined by eqs.( p.llD as 

a = e'^'ustD^i^p), a = e'^'u^.tDgip), 
are equal to each other 

a = a= ^^^'^jJ^^^K ^lI^^ , (5.7) 

where the subscripts denote partial differentiations. Hence the necessary 
condition for a solution to be conformally invariant is the equality 

a = a ( ^ A = A). (5.8) 

The converse statement gives the criterion for a solution to be conformally 
non- invariant. 

Corollary 1 The sufficient condition for a solution of the heavenly equation 
to be conformally non-invariant is that the following inequality should be 
satisfied 

a^a (5.9) 

(or equivalently A 7^ A). 

Concerning the practical use of this statement we must remark that even 
if the inequahty ( ^.9|) is satisfied for a solution of the resolving system it 
could become the equahty ( |5.8| ) on the corresponding solution of the heavenly 
equation. Nevertheless the above criterion is useful, meaning that we should 
avoid solutions of the resolving equations satisfying eq.(|5.8|) in order not to 
end up with conformally invariant solutions. 
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6 Particular solutions of the resolving 
system 

Here we show that the commutator representation of the resolving system can 
prompt Ansatzes, leading to particular solutions of the resolving equations. 
Attempts to solve the commutation relations by imposing relations between 
the operators of invariant differentiation lead to invariant solutions of the 
heavenly equation. This is the case with the Ansatz Y = Y. Then the 
expressions ([4.11 ) for Y, Y imply A = A, so the condition (|5.9| ) of Corollary 1 
is not satisfied. Hence we obtain a conformally invariant solution of the 
heavenly equation. 

Another possible simplifying Ansatz is that the operators Y and Y com- 
mute and we have 

T = -utp [r,F]=o, (6.1) 

but Y ^Y. 

Before solving the resolving system with the Ansatz ( |6.1| ) we keep in 
mind that F ^ 0. Indeed the case F = is singular for the derivation of the 
resolving equations and should be treated separately. We shall consider first 
the case F = and show that it leads to invariant solutions of the heavenly 
equation. 

Putting F = in equation ( |4.1|) we obtain 

Uzt = 0, u,t = (6.2) 
and hence we have the separation 

u = a{t) + f3{z,z). (6.3) 
Substituting this expression for u into the heavenly equation ( |2.1| ) we obtain 
{a"{t) + {a'it)f) = ne-^^^''^ f3,,{z, z) = 21, (6.4) 

where / = / is a separation constant and the primes denote derivatives in t. 
Integrating the equation for a we obtain 

a{t) = In (it^ + Cit + C2) (6.5) 

where Ci,C2 are arbitrary real constants. The equation for (3{z,z) 

(3,-, = 2K/e^ (6.6) 
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is the Liouville equation if k = 1 and the "pseudo-Liouville' equation for 
K = —1. Its general solution is 

f3{z, z) = In a'{z) + In a'{z) - 2 In {a{z) + a{z)) - In / (6.7) 

if K = 1 and 

(3{z, z) = In a'{z) + In a{z) - 2 In {a{z)a{z) + 1) - In / (6.8) 

ii K = —1. Here a{z) is an arbitrary holomorphic function and the primes de- 
note derivatives. Thus, the corresponding solutions of the heavenly equation 
are given by the equation ( |6.3| ) with a(t) determined by the formula ( |6.5| ) 



and f3{z,z) determined by the formula ([6171) or (|6l^ ). 

To obtain the simplest representative of the orbit of solutions we apply 
simplifying symmetry transformations: the conformal transformation 

a{z) z, d{z) (— > z, 

the suitable time translation and the dilation of time accompanied by a shift 
of u 

M I— >• M + In /, t 



t 



The resulting solutions become 

u = \n(t^ + C) -2\n{z + z) if k = 1, (6.9) 



u 



\Yi(f + C) -2\n{zz + l) if /€ = -!, (6.10) 

where C is an arbitrary real constant. 

To perform a check of invariance of the solutions ( p.9|) and (|6.10|) , we 
substitute them into the criterion of invariance ( ^.3|) and make a splitting in 
t. Then we obtain a = and if C 7^ 0, then also /3 = 0. For C = the 
constant j3 can be arbitrary. We also obtain a differential equation for a{z) 
and a{z) 

a'{z) + -a'{z)=2^^^^^^^^ for k = 1 (6.11) 

z + z 

and 

i( \ , -'(-\ ^za{z) + za{z) 
a[z) + a[z) = 2 —^—^ for k = — 1, (6.12) 
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with the trivial solutions 

a = i, a = —i for k = 1; a = iz, a = —iz for k = —1. (6.13) 

Thus, we have proved that there exist a,/3 and a{z),a{z) such that the cri- 
terion ( p.3|) of invariance of solutions is satisfied for our solutions ( |6.9| ) and 
( 6.101 ). Hence the case F = corresponds to invariant solutions. 

In the following we assume that F 7^ and consider the resolving equa- 
tions with the Ansatz (|6.1|) . Equations ( [4. 61) and (|4.1CI| ) become respectively 

F(A) = F(A) and Y{X) + F(A) = 4k 

and hence 

F(A) = 2k, F(A) = 2k. (6.14) 

Next we consider the compatibility condition of the system of equations ( [4.8|) , 
(|4.9| ) and the first equation in ( |6.14| ). Because of the formula (|6.1|) the first 
equation becomes 

5{\) =ut\- - p. (6.15) 

Then, using cross-differentiation of the invariant operators 5 and Y , their 
commutator ( 3.14 ) and eq.( 6.14 ), we obtain a very simple result 

A + A = 2KMi. (6.16) 



Solving this equation with respect to A, substituting in equation ( |4.9| ) and 
using the equation ( |6.15|) to express 5(A), we obtain a quadratic equation for 
A 

A^ - 2KUt\ + 2kp = 0, 

with the solution 



X = KUt + iy2Kp — uf, (6-17) 

where we have chosen the -|- sign before the square root. Equation (|6.16|) 
gives the result 



X = KUt — iy 2kp — Ut (6.18) 



which is complex conjugate to ( |6.17|) provided the condition 

2kp-m^>o (6.19) 

is satisfied. The obtained expressions for A and A satisfy the equations ( p.l4| ), 
( 6.15 ) and its complex conjugate, and hence all the resolving equations apart 
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from the equation ( [4.4| ). We rewrite this last equation for the new unknown 
/ = In F as 

ft + {np - uDfu, - Utpfp = -3ut (6.20) 

and solve it by the method of characteristics obtaining the general solution 
of the equation ( [4.4|) 

F = pV(e, 0) where ^ = e = t-- (ut + ^2Kp-u'i^ , 

^ ^ (6.21) 

where 99 is an arbitrary real differentiable function. 

Finally we sum up our results for the particular solution of the resolving 
system which follows from our Ansatz (|6.1|) 

F = p^ip{^,9), T = —utp, A = KUt + i\J 2kp — Ut, A = nut — i\J 2Kp — 

(6.22) 

This will be used in the next section for obtaining the corresponding solution 
of the heavenly equation. Since A 7^ A the condition ( ^.9|) of Corollary 1 for 
non-invariance of this solution is satisfied. 



7 Invariant integration and non-invariant 
solution of the heavenly equation 

In this section we reconstruct the solution of the heavenly equation starting 
from the particular solution (|6.22|) of the resolving system. We demonstrate 



here the procedure of invariant integration which amounts to the transfor- 
mation of equations to the form of the exact invariant derivative. Then we 
drop the operator of invariant differentiation adding the term playing the 
role of the integration constant which is an arbitrary element of the kernel 
of this operator. 

We start from our Ansatz (|6.1|) using the definitions t = 6{p) and 6 = Dt 



Dt{\np) = Dt{-u). (7.1) 
We integrate this equation in the form 

Inp = -M + \ia-f^g{z,z), 
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where the last term is a function to be determined. Solving this equation 
with respect to p and using the definition of p we obtain 



p = e Uzz = e 'yzz{z,z) 
and hence Uzz = lzz{z, z). This implies the following form of the solution 

u{z, t) = 7(2;, z) + a{z, t) + d{z, t), (7.2) 

where 7, a and a are arbitrary smooth functions of two variables. After the 
substitution of this expression into the heavenly equation it becomes 



^a(z,t)+a{z,t) 



au{z, t) + au{z, t) + {at{z, t) + at{z, t))" 



(7.3) 

Next we rewrite the formulae ( |6.22 ) for A and A in the form of exact 
invariant derivatives 



Y{J2kp — u1 — inut) = 0, Y{J2kp — u1 + iKUt 



0. 



(7.4) 



On account of the definitions (|3.13| ), the operators Y and Y can be written 
as 



and the equations ( [7.4|) become 



'2np — Ui + iK,Ut)z = 0, 
They are integrated in the form 



2K,p — u1 — tKUt) 



2kp — + inut = "ipiz, t), v 2fi;p — ~ ihiUt = ip{z, t) 



(7.5) 



where ip, ip are arbitrary smooth functions. Taking the difference of two 
equations ([7.5|) we obtain 



IK 



Ut 



— [i){z, t) - -^{z, t)] = at{z, t) + at{z, t), 



where the last equality follows from the expression ( |7^ ) for u. Separation of 
z, z in the last equality leads to 



IK 



2^ 

atiz,t) - —^p{z,t) 



x'it) = -x'it), 
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where x'{t) is the separation 'constant" and the prime denotes the derivative. 
Solving these equations with respect to ip, and substituting the results into 
the equations ( |7.5| ) we solve them with respect to the square root with the 
result 

\j2Kp - uf = in [at{z, t) - at{z, t) - 2x'(t)] • 

Solving this equation with respect to up and multiplying the result by 6°"^°" 
we obtain 

ne-^^'-'hz-z{z,z) (7.6) 
= 2e"(^'*)+^(^"'*) [at{z,t)at{z,t) + - at{z,t)) - x'\t) . 

Using this equation in the right-hand side of the heavenly equation in the 
form (|773|) and separating z, z we obtain two complex conjugate equations 

au{z, t) = -a\{z, t) + 2x'{t)at{z, t) - x'\t) + fi{t), (7.7) 
au{z, t) = -al{z, t) - 2x'{t)at{z, t) - x'\t) - ^{t), (7.8) 

where p(t) = —pit) is the separation "constant". We substitute these ex- 
pressions for att and au into the transformed heavenly equation (|7.3|) to 
obtain 



atiz, t)at{z, t) + x'{t){at{z, t) - atiz, t) - x'^(t)) 



^e-^^'~'hAzrz). (7.9) 

Next we take the total derivative Dt of this equation and substitute again 
the second derivatives au and otu from the equations (|7. 7|) and (|7.8| ). The 
result is unexpectedly simple 

(X" - /x) («t - - 2x0 = 0. (7.10) 

This equation implies that 

p{t) = x"{t), (7.11) 
since the complementary assumption 

at - at - 2x' = 

leads again to the equation (|7.11|) . Indeed, the last equation allows a sepa- 
ration of z, z and, being integrated, gives a, a 

a{z, t) = x{t) + iy{t) + uj{z), a{z, t) = -x{t) + u{t) + Q{z). 
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Substituting these expressions into the equations ( [7. 71 ) and ( |7.8| ) and com- 
paring the results we discover again the equation (|7.11|) . 

With this restriction the equations ( |7.7] ) and ( |7.8| ) are simphfied and 
integrated to give 

a{z, t) = In (t + b{z)) + x{t) + uj{z), a{z, t) = In (t + a{z)) - + co{z), 

(7.12) 

where b{z) and uj{z) are arbitrary holomorphic functions and we have reserved 
the notation a{z) only for the generators of the conformal vector field Xa in 
eq.(|3). 

Now define a new function of z, z 

T{z,z) =-f{z,z)+u{z)+iu{z), (7.13) 

so that the form ( [7.2| ) of the solution becomes 

u{z, z, t) = In (t + b{z)) + In (t + b{z)) + r{z, z). (7.14) 

Substituting the expressions (|7.12|) for a, a into the transformed heavenly 
equation ( |7.3| ) we obtain the equation for the only unknown function T{z,z) 
in the solution ( [7.14| ) 

T,, = 2Ke^. (7.15) 
If K = 1 this is the Liouville equation with the general solution 

V{z,z) =\nc{z)+\nc\z)-2\Yi{c{z)+c{z)) (7.16) 

where c{z) is an arbitrary holomorphic function. If k = — 1 we call the 
equation ( |7.15| ) 'pseudo-Liouville' equation and its general solution is 

T{z,z) =lnc'(z)+lnc'(z)-21n(c(^)c(z) + l). (7.17) 

Finally, substituting these expressions for r(z, z) into the equation ( |7. 14| ) we 
obtain the solutions of the heavenly equation ( |2.1j ) for the two choices of the 
sign K = +1 and k = —1. 

1. The solution for k = 1: 

u{z, z, t) = In (t + b{z)) + In (t + b{z)) 
+ \nc\z)+\nc\z)-2\n{c{z)+c{z)). (7.18) 
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2. The solution for k 



-1 (see also 



u{z, z, t) = ln{t + b{z)) + In (t + b{z)) 
+ In c{z) + lnc'{z) - 2 In {c{z)c{z) + 1). 



(7.19) 



Here b{z) and c{z) are arbitrary holomorphic functions. 

To avoid "false generality" it is sufficient to choose the simplest represen- 
tative of the obtained orbits of solutions applying the conformal symmetry 
transformation c{z) = z, c{z) = z with the following results. 

1. The solution for k, = 1: 

u{z, z, t) = ln{t + b{z)) + In (t + b{z)) - 2 In + z). (7.20) 



2. The solution for k = —1: 



u{z, z, t) = In (t + b{z)) + In (t + b{z)) - 2 In {zz + 1). (7.21) 



Here b{z) is still an arbitrary holomorphic function. 

Up to now we solved completely only the Ansatz ( |6.1| ) defining r, but 
we did not check the automorphic equation ( [4.1|) and the auxiliary equations 
(|4.5|), by using the particular solution ( |6.22D of the resolving system. Hence, 
though we obtained the correct solutions ( [7.20| ) and (|7.21| ) of the heavenly 
equation ( p.l|) , we have not made a complete foliation of these solutions into 
separate orbits. 

To do this, first we remark that due to the discrete symmetry of the 
heavenly equation ( |2.1| ) and of our solutions with respect to the permutation 
z ^ z,we can define the holomorphic function b{z) as satisfying the condition 



lmb{z) > for K = 1 and lmb{z) < for k = —1. 



(7.22) 



Then we check that the automorphic equation (|4.1| ) coincides with the aux- 
iliary equations (^4.5|) and becomes 



{z + zYb'{z)b'{z) = 8^{^,i 



for K = 1 



(7.23) 



and 



{zz + l)%'{z)b'{ 



"z) 



-8^(e, 



for 



K 



■1. 



(7.24) 
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Using the solutions ( [7.2(J| ) and ( |7.21| ) in the definitions ( |6.21| ) of the char- 
acteristic variables ^ and 9^ we discover that they depend only on h and 6, 

^_(^_ ^_(^_^ ^^^^^ 

Hence, defining the new arbitrary function $(6, h) = ip{C,, 6), the automorphic 
equations ( [7 .231) and (|7.24| ) become 

{z + zfb'{z)b'{z) = 8$(6, b) for k = 1 (7.26) 

and 

{zz + l)%'{z)b'{z) = -8^{b, b) for k = -1. (7.27) 

Sufficient conditions for solving these functional-differential equations are 
given by the following choices of $(6, b) 

'fib) + f{b)Y 

<^{b,b) = ^ V/, for K=l (7.28) 



and 



$(6, b) 



if'{b)r{b) 



f{b)f{b) + 1 



for K, 



(7.29) 



8/'(6)/'(6) 

where f{b) is an arbitrary holomorphic function. Then the automorphic 
equations become 



{in [fib) + m]}^_^ = [\niz + z)l 



for K = l (7.30) 



and 



{in [f{b)m + l]}^^=[\n{zz + l)l. 
Their general solutions are 

/(6) + /(6) = u;(z)^(^)(z + ^) 

and 

f{b)f{b) + 1 = w{z)w{z){zz + 1) 

where w{z) is an arbitrary holomorphic function. The formulae ( [7. 281 ) and 
( |7.29| ) for $(6,6) become 



for K = -1. (7.31) 



for K = 1 (7.32) 



for K = -1, (7.33) 



$(6,6) 



w'^{z)w'^{z){z + zY 
8f(6)/'(6) 
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for K = 1 



and 

= — smm ^ = -^' 

If we plug these formulae into the automorphic equations (|7.26|) and ( [7.27|) , 
then both automorphic equations coincide and become 

n.m.)^^. (7.34, 

This equation admits a separation of variables, leading to the ODEs 

The obvious choice of the functions w{z) and w{z) 

w{z) = 1 w{z) = 1 (7.36) 
simplifies the ODEs ( fTSSD to 

[f{b)l = 1, = 1, (7.37) 

with the solution 

mz)] = z, ]\h{z)\ = z (7.38) 

meaning that 6(2;) is the inverse function for /(&): h = f^^. The equations 
( |7.32| ) and ( |7.33| ) are obviously satisfied by the solution ( [7.38| ) with our choice 
(17361 ) oiw{z),w{z). 

Thus, any particular function b(z) can be obtained for an appropriate 
choice of f{b) as its inverse function. This fixes the function $(6, b) according 
to the formulae ( |7.28| ) or ( |7.29| ), the function (p{^, 6) = $(6, b) and the right- 
hand side F of the automorphic equation (|4.1|) determined by the formulae 
( 6.22 ). Hence any particular choice of the function b{z) in our solutions 
( |7.2CI| ), ( |7.21| ) means a corresponding choice of the particular orbit in the 
solution space of the heavenly equation. 



8 Check of non- invar iance of the solutions 



In this section we prove that our solutions ( [7.20|) and ( |7.21| ) of the heavenly 
equation (|2.1| ) are non- invariant, with respect to its symmetry group gener- 
ated by the vector fields in ( |2.2| ), for generic functions a{z)^ except for some 
particular classes listed below in the theorems summarizing the results. 
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For the check of non-invariance we substitute our solutions ( |7.20| ) and 
(|7.21| ) into the invariance criterion (|5.3|) . The resulting equation is quadratic 
in t and it implies the vanishing of the coefficients oi t^,t and 

We consider first the case k = 1. The term with gives again the 
equation (|6.11|) of the Section 6. However, now we need the general solution 
of this equation. 

We assume in the generic case that a'{z) + a'{z) ^ 0, otherwise the 
equation ( |6.11| ) implies a = —a = constant and this case should be treated 
separately. Then we rewrite the equation (|6.11|) in the form 



a{z) + a{z) 
a'(z) + a'(z) 



z + z 



a + a 
a' + d' 



0. 



(8.1) 



In order to consider the generic case we postulate a"a" ^ 0, then the last 
equation can be easily manipulated, obtaining the solution 



a\z] 



Ci(^ + A)2 + C2, 



a\z] = — 



Ciiz-\f + C2 



(8.2) 



where Ci ^ 0, C2 and A are arbitrary purely imaginary constants. 

Now we consider the term without t in the criterion of invariance using 
our result (^.2|) which gives the equation with the separated variables z, z 



a 



h{z] 



+ 



C^{z + \f + C2 



h'{z) 



h{z) 



(8.3) 



a 

W) 



Ci{z-\f + C2 



m 



/3 



where is a separation constant. Comparing these equations with the equa- 
tion obtained from the term with t in the criterion of invariance we conclude 
that they coincide if and only if the condition 



^^ {b{z) - b{-z) 







is satisfied. This implies /x = 0, since otherwise we have b = b = constant 
and our solution is obviously invariant depending only on two variables t and 
z + z. Hence the equations (^.31) become 



(8.4) 
(8.5) 



Ci{z + + C2] b'iz) - (3biz) = -a, 
Ci{z - A)2 + C2I b'{z) + f3b{z) = a. 



23 



Consider now the case C2 7^ 0, /? 7^ and introduce the notation 



Integrating the ODEs ( ^.4| ) and ( ^.5| ) we fix the function b{z) in our solution 
of the heavenly equation which corresponds to the invariant solution in the 
considered case 

h{z) = C ^ + -, b{z) = C - — ^ + - 8.7 

where C, C are integration constants. 

In a similar way we treat other possible cases. We sum up the results for 
the case of k = 1 in the following statement. 

Theorem 2 The function 

u = hi{t + h{z)) + In (t + h{z)) -2\yi{z + z) (8.8) 

is a solution of the heavenly equation for k = +1 for an arbitrary 

holomorphic function b{z). This solution is a non-invariant solution of this 
equation iff the function b{z) does not coincide with any of the following 
choices: 

1. 

\z + X + u ) (3 

where a and (3 are arbitrary real constants, (3 0, u and 7 are defined 
by the formulae 6| j and A, Ci, C2 are complex constants which satisfy 
the conditions 

X = -X, Ci = -Ci, C2 = -C2, Ci^O, C2^0. 

In this case the solution is invariant with respect to the symmetry gen- 
erator 

X = adt + (3 {tdt + 2du) + Ci [{z + Xfd, -{z- Xfd, - 2{z - z)d^ 
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2. 



biz) 



a f z + X + 
2V-C1C2 U + A-i/^ 



+ C if /3 = 0, C2 ^ 0; 



t/ie solution is invariant with respect to the previous symmetry generator 
X with /3 = 0. 



3. 



b{z) — C exp 



/3 



Ci(^ + A) 



+ 



a 



if C2 = 0, /3 7^ 0; 



the solution is invariant with respect to the symmetry generator X from 
the case 1 with C2 = 0. 



b{z)^ 



a 



Ci(z + A) 



+ C if /? = and C2 = 0; 



the solution is invariant with respect to the symmetry generator X from 
the case 1 with (3 — and C2 — 0. 



5. 



b{z)^C{C,z + C2f'^' + '^ if Ci^O, /3 7^0; 

the solution is invariant with respect to the symmetry generator 
X = adt + p {tdt + 2d^) + Ci {zd, + zd, - 2aj + C2 {d, - d,) . 



6. 



7. 



13 (y 

b{z) = Ce-2' + - if Ci = 0, /5 7^ 0; 

the solution is invariant with respect to the symmetry generator X from 
the case 5 with Ci — 0. 



h{z) = -^z + C if Ci = 0, /3 = 0; 

the solution is invariant with respect to the symmetry generator X from 
the case 5 with Ci — and P — 0. 
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8. 

b{z) = b = constant if Ci = a = (3 = 0, C2 0; 
the solution is invariant with respect to the symmetry generator 

X = d,- d,. 

If b = a/P, then this solution is also invariant with respect to the 
generator 

X = adt + [3 {tdt + 2du) . 

Now we consider the 1 and substitute the solution (|7.21|) of 

the heavenly equation into the criterion of invariance ( |5.3| ). Then the 
resulting equation is again quadratic in t and the term with gives us again 
the equation ( |6.12| ), for which we need now the general solution. First we 
rewrite it in the form 

, , a' + a' ^ f a a\ 
a +a! + — = 2 - + - . 

zz \z zj 

Differentiating this equation with respect to z and z we obtain an equation 
which admits separation of z, z in the form 

zal'^z) - a\z) = - [za"{z) - a'{z)] = A = -A (8.9) 

where A is a separation constant. Integrating these ODEs we obtain 

a{z) = Ciz"^ - Xz + C2, a{z) = Ciz"^ + Xz + C2 

where Ci, C2 are integration constants. Substituting these solutions into the 
equation (|6.12|) we see that it is identically satisfied if and only if 
Ci = C2 <^==^ 6*2 = Ci, so that finally we have the solution of the equation 
following from the term with t^ 

a{z) = Ciz^ - Xz + C2, a{z) = C2z'^ + Xz + Ci. (8.10) 

Next we consider the term without t in the criterion of invariance using 
our result ( ^.101) which gives the equation with the separated variables z, z 



+ {C,z^-Xz + C2)"-^^P 



b{z) V ^ 'V 

>'(^\ 1 

/i = -/i (8.11) 



" +(C2z' + Xz + cAI^-P 



b{z) ' V ^ ' V 
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where /i is a separation constant. Comparing these equations with the equa- 
tion obtained from the term with t in the criterion of invariance we conclude 
that they coincide if and only if the condition 



/i 



(b{z) - b{z)) = 



is satisfied. This implies /i = for the same reason as in the case k = 1. 
Hence the equations become 

[Ciz^ -Xz + C2) h\z) - Pb{z) = -a, (8.12) 

[C2Z^ + Xz + Ci) b'{z) - (3b{z) = -a. (8.13) 

Consider now the case Ci 7^ 0. Introduce the new constants A = — A/(2Ci) 
and C2 = C2 — A^/ (4Ci). Then the first equation takes the form 

Ci{z + ~Xf + 02] b'{z) - (3b{z) = -a (8.14) 



coinciding with the ODE (|8.4| ) in the case k = 1. Therefore we can use its 
solutions with an appropriate change of notation. Other possible cases are 
treated in a similar way. Therefore we can transfer the results of Theorem 
2 to the case k = —1 with an appropriate change of notation and sum them 
up in the following statement. 

Theorem 3 The function 

M = In (t + b{z)) + In (t + b{z)) - 2 In {zz + 1) (8.15) 

is a solution of the heavenly equation ( |<^. i[ ) for k, = —1 for an arbitrary 
holomorphic function b{z) . This solution is a non-invariant solution of this 
equation iff the function b{z) does not coincide with any of the 8 forms given 
in Theorem 2 with the change of notation 

A ~ A^ 

X ^ X = -— , C2 6*2 = Co 



2Ci 4Ci 



C2 - 13 

M-^, 7-7 = 



C1C2 



in the cases 1, 2, 3, 4 and Ci ^ —X in the case 5. Those 8 choices ofb{z) give 
invariant solutions with respect to the corresponding symmetry generators of 
Theorem 2 with the same change of notation. 
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9 Conclusions and outlook 



The title of this article, or rather of the research direction that it represents, 
could have been "Invariant methods for obtaining non-invariant solutions 
of partial differential equations". The main result is that we are proposing 
an alternative tool for obtaining particular solutions of non-linear partial 
differential equations with infinite dimensional symmetry algebras. As stated 
in the Introduction, the idea of the method is more than a hundred years 
old We have turned it into a usable tool by adding new elements. 

These are: 

1. The systematic use of invariant cross-differentiation involving the oper- 
ators of invariant differentiation and their commutator algebra for the 
derivation of the resolving equations and for obtaining their particular 
solutions. 

2. The presentation of the resolving system as a Lie algebra of the opera- 
tors of invariant differentiation (over the field of differential invariants 
of the symmetry group) |jTO|. 



3. The concept of invariant integration applied to the automorphic system. 



Let us use the heavenly equation (|2.1|) to compare different methods of 
obtaining exact analytical solutions of a partial differential equation, pro- 
vided or at least suggested by symmetry analysis. In all of them the studied 
equation is embedded into a larger system of equations, to be solved simul- 
taneously. 

The most standard method is that of invariant solutions [^,|l^,[l6|. One 



first finds the symmetry algebra realized by vector fields of the form 

x = Tdt + id, + ld, + ct>du (9.1) 

where r, ^ and are functions of t, z, z and u. Once this algebra is found 
(i.e. the algebra ( |2.2| ) for the heavenly equation) one classifies its subalgebras 
into conjugacy classes and then adds one, or more, first order linear equations 
of the type 

rut + iu, + - = (9.2) 

to the studied equation. These equations are solved, their solution is substi- 
tuted into the original equation. This again is solved and we obtain solutions 
invariant under the chosen subgroup. 



28 



Further methods are the Bluman and Cole "non-classical method" IT? 



the Clarkson-Kruskal "direct method" and that of "conditional symme- 
tries" ||19[ (see for a review). These methods, basically all equivalent, 
amount to the fact that a first order equation of the type (|9^) is added to 
the studied equation, without the requirement that t,^,C, and (j) define an 
element of the symmetry algebra. 



Finally, we have the group foliation method |jTO[ used and further devel- 



oped in this article. Let us review the essential steps, performed above. 
1. Find the total symmetry algebra (|2.2|). 



2. Find all differential invariants of order up to of its infinite dimen- 
sional subalgebra which is Lie algebra of the conformal group. The 
number must be larger or equal to the order of the equation and 
must satisfy the requirement that there should be functionally 
independent invariants with 

i^N>p + q (9.3) 

where p and q are the number of independent and dependent variables, 
respectively. In our case we have p = 3, q' = 1, N = 2, = 5. The 
actual invariants are given in the equation 



Choose p invariants as new independent variables and require that the 
remaining invariants be functions of the chosen ones. This provides us 
with the automorphic system that also contains the considered equa- 
tion, expressed in terms of the invariants. In our case the automorphic 
system consists of the equation ( p.8|) (the heavenly equation) and the 
equation (|4.1j ) (or equivalently (|4.2| )). 



Find the "resolving equations". This is a set of compatibility condi- 
tions between the studied equation and those that we have added to 
obtain the automorphic system. In our case we require compatibility 
between the equations ( |2.8| ) and (^?T|), i.e. determine the restrictions 
on the function F{t,Ut,p). We have shown that this can be done in 
an explicitly invariant manner by using the operators of invariant dif- 
ferentiation, in our case S,Y and Y of the equations (|3^ ) and ( p.l3|) . 



The resolving system in our case consists of the equations (|4.4| ), ( [4.6|) , 
( [4.9| ) and ( [4.10| ). As stated by the fundamental Theorem 1, this 



resolving system is best written as a system of commutator relations 
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for the operators of invariant differentiation projected on the solution 
manifold of the heavenly equation in the space of differential invariants, 
together with the Jacobi relations for these operators. 

5. Solve the resolving system and the automorphic one. This provides 
solutions of the original equation. 

The last step, step 5 is the most difficult one. If it can be carried out 
completely, we obtain "all" solutions, both invariant and non-invariant ones. 
In general, such a situation is too good to be true. In particular, for the heav- 
enly equation we were not able to solve the system ( |3.14|) , ([4.1 3|) in general. 
Instead, we made various simplifying assumptions. The most obvious ones, 
like Y = Y or F = 0, lead to invariant solutions. These we already know, 
or can obtain by much simpler standard methods. The assumption, or re- 
striction, that leads to non-invariant solutions was [Y, Y] = 0. The solutions 
obtained are ( [7.18|) and (|7.19|) , for k = 1 and k = —1, respectively. Each 



solution involves two arbitrary holomorphic functions. One of them, b{z) is 
fundamental. The other is induced by a conformal transformation and can 
be transformed away {i.e. set equal to e.g. c{z) = z). In Section ^ we show 
that the solutions are, in general, not invariant under any subgroup of the 
symmetry group. They reduce to invariant ones only for very special choices 
of the function h{z), specified in Theorems 2 and 3. 

It would be interesting to relate the concepts of this article to that of inte- 
grability for non-linear partial differential equations. "Integrability" means 
that the considered equation is viewed as an integrability condition for a 
Lax pair, a pair of linear operators [^T],|2^. Here we can view the equations 
(|3.12| ) as a set of relations between a triplet of linear operators, subject to a 



non-linear constraint ( [4.13| ) . 
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